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HIGHLIGHTS 


•  Thermoelectric  modules  with  various  leg  geometries  were  modeled  and  analyzed  for  two  different  temperature  gradients. 

•  Disparities  in  temperature  distributions,  power  outputs,  and  conversion  efficiencies  were  limited. 

•  Magnitudes  and  distributions  of  thermal  stresses  in  the  legs  were  significantly  affected  by  changing  leg  geometries. 

•  Thermal  stress  levels  were  smaller  in  the  cylindrical  and  trapezoidal  prism  legs  compare  to  the  rectangular  prism  legs. 
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This  study  aims  to  investigate  possible  effect  of  various  thermoelectric  leg  geometries  on  thermo¬ 
mechanical  and  power  generation  performances  of  thermoelectric  devices.  For  this  purpose,  thermo¬ 
electric  modules  with  various  leg  geometries  were  modeled  and  finite-element  analyses  for  two  different 
temperature  gradients  were  carried  out  using  ANSYS.  Temperature  distributions,  power  outputs,  con¬ 
version  efficiencies,  thermal  stresses  in  the  legs  were  evaluated  for  each  model.  Significant  differences  in 
magnitudes  and  distributions  of  the  thermal  stresses  in  the  legs  occurred  due  to  changing  leg  geome¬ 
tries.  Thermal  stresses  in  the  rectangular  prism  and  the  cylindrical  legs  were  49.9  MPa  and  43.3  MPa 
respectively  for  the  temperature  gradient  of  100  °C. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Environmental  challenges,  such  as  global  warming,  growing 
demand  on  energy,  and  diminishing  oil  sources  have  accelerated 
research  on  alternative  energy  conversion  methods.  Meanwhile, 
significant  amount  of  energy  is  wasted  during  energy  conversion 
processes  due  to  technological  limitations.  Particularly  in  vehicles, 
more  than  50%  of  the  fuel  energy  [1  ]  dissipates  to  the  environment 
through  exhaust  and  cooling  systems  as  heat.  Recovering  some 
amount  of  wasted  heat  and  converting  it  into  useful  energy  such  as 
electricity  could  increase  overall  efficiency  of  energy  conversion 
systems  and  reduce  demand  on  fossil  fuels  and  natural  resources. 
Thermoelectric  power  generation  is  an  appropriate  and  promising 
method  to  convert  wasted  heat  into  electrical  energy.  Imposed 
temperature  gradient  through  a  thermoelectric  device  (module) 
produces  electrical  potential  through  the  Seebeck  effect,  discovered 
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by  Thomas  J.  Seebeck.  Thermoelectric  devices  can  also  be  used  for 
cooling/heating  purposes  utilizing  Peltier  effect. 

Based  on  the  application,  thermoelectric  devices  can  be  used  for 
cooling,  heating,  power  generation,  and  sensing  2].  The  benefits 
include  being  solid  state  (no  moving  parts),  silent,  maintenance 
free  and  light  weight.  These  characteristics  attract  attention  to 
thermoelectric  devices  from  industries  including  defense,  aero¬ 
space,  electronics,  and  automotive.  Temperature  controlled  car 
seats,  electronic  device  cooling,  thermal  management,  waste  heat 
recovery  from  hot  chimneys,  thermoelectric  cogeneration  systems 
and  energy  harvesting  from  human  body  heat  are  some  of  the 
existing  and  potential  future  applications  of  thermoelectric  systems 
[1-10].  Thermoelectric  generators  (TEGs)  have  the  capability  to 
convert  wasted  exhaust  heat  into  useful  electric  power  in  vehicles 
[11].  Thus,  automotive  companies  have  begun  developing  pro¬ 
totypes.  Nonetheless,  thermoelectric  devices  have  not  been  broadly 
used  particularly  for  power  generation  due  to  some  of  their  de¬ 
ficiencies.  Low  conversion  efficiencies,  usually  around  3—10%  [1], 
due  to  the  small  non-dimensional  figure  of  merit  ZT  value,  is  the 
major  drawback  of  thermoelectric  devices.  New  thermoelectric 
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materials  (with  values  of  ZT~1  and  more)  including  thin-film 
superlattices,  advanced  bulk  semiconductor  alloys,  and  crystal 
structures  with  rattlers  have  been  developed  recently  through 
intensive  research  efforts  focused  on  enhancing  efficiencies  [12]. 
Contrary  to  the  material  level  challenges,  research  on  device-level 
challenges  is  small.  This  resulted  as  an  imbalance  between  mate¬ 
rial  and  device  level  developments  of  thermoelectric  devices. 

Thermoelectric  devices  are  typically  consists  of  electrically  in 
series  and  thermally  in  parallel  connected  p-  and  n-type  thermo¬ 
electric  legs,  commonly  fabricated  with  BiTe-,  PbTe-  or  SiGe-based 
semiconductors.  Conductor  shunts,  cold/hot  substrates  and  solder 
joints  are  the  other  parts  of  a  thermoelectric  device.  In  either  power 
generation  or  cooling  applications,  they  are  subjected  to  a  tem¬ 
perature  gradient.  This  causes  thermal  stresses  in  the  devices  due 
to  differential  thermal  expansions  and  mismatching  of  thermal 
expansions  of  the  bonded  components  [13  .  Dislocations  and  cracks 
might  arise  due  to  the  thermal  stresses.  Further,  fatigue  fractures 
may  occur  and  finally  the  device  may  fail  due  to  repeated  thermal 
stresses  in  cyclic  heating  applications  [14,15].  Potential  failure 
mechanisms  of  thermoelectric  devices  were  evaluated  by  Choi  et  al. 
[15].  More  seriously,  failure  of  even  only  one  thermoelectric 
element  may  cause  the  entire  device  to  fail  since  all  the  elements 
are  connected  electrically  in  series.  Consequently,  thermal  stresses 
are  critically  important  for  reliability;  operation  life-time  and  ex¬ 
pected  performance  of  thermoelectric  devices.  Therefore,  it  is 
crucial  to  investigate  developed  thermal  stresses  in  thermoelectric 
devices  and  the  factors  that  influence  them. 

Researchers  carried  out  numerical  analyses  to  determine  power 
generation  and  thermo-mechanical  characteristics  of  thermoelec¬ 
tric  devices  and  effect  of  geometric  factors.  Russel  and  co-workers 
[8,9],  numerically  and  experimentally  studied  performance  char¬ 
acterization  of  a  thermoelectric  cooling-based  (TEC)  thermal 
management  system  while  considering  parameters  include  geo¬ 
metric  factor  of  the  thermo-elements.  Jang  and  co-workers  [5,6], 
studied  power  generation  performance  of  thermoelectric  modules 
using  numerical  methods  to  find  optimum  performance  parame¬ 
ters.  Li  and  co-workers  [13]  showed  that  the  thickness  of  the  copper 
pad  does  not  have  an  effect  on  the  power  output  of  the  thermo¬ 
electric  device  under  defined  conditions.  Jiang  and  co-workers  [16] 
studied  effect  of  thermoelectric  leg  size  and  variable  physical  pa¬ 
rameters  on  power  generation  and  conversion  efficiencies  of  ther¬ 
moelectric  modules.  Furthermore,  Nguyen  and  Pochiraju  [17] 
studied  power  generation  performance  of  thermoelectric  modules 
under  transient  temperature  gradient  conditions. 

Studies  on  thermo-mechanical  performance  revealed  that  the 
maximum  thermal  stresses  are  developed  at  the  interface  between 
the  hot  side  conductor  and  the  legs  of  BiTe-based  modules 
[14.18,19].  The  corners  of  the  legs  are  more  critical  through  the 
concentrated  stress  levels.  Researchers  showed  that  varying  leg 
lengths  and  widths  changes  the  thermal  stresses  developed  in  the 
modules  [19-21].  Nakatani  and  co-workers  [22  observed  that 
increasing  conductor  thickness  lowers  the  stress  caused  by 
macroscopic  deformation  despite  increases  in  the  strain  caused  by 
CTE  mismatches.  The  influence  of  various  trapezoidal  prism  leg 
geometry  configurations  on  thermal  stress,  power  generation,  and 
the  device  efficiency  of  the  thermoelectric  modules  was  also 
studied  [23,24].  Al-Merbati  and  co-workers  [23]  showed  that  lower 
thermal  stresses  can  be  acquired  by  changing  leg  geometry  and  the 
consequence  of  this  can  enhance  operation  life.  Investigation  led  by 
researchers  showed  distribution  of  thermal  stresses  in  segmented 
thermoelectric  legs  25-27].  Li  and  co-workers  [25]  observed 
maximum  tensile  stress  of  180  MPa  due  to  mismatching  of  thermo¬ 
mechanical  properties.  Consequently,  thermal  stresses  are  critically 
important  for  durability,  operation  life-time,  and  expected  perfor¬ 
mance  of  thermoelectric  devices.  We  agree  with  several 


researchers  [28,29],  that  thermoelectrics  needs  more  device-level 
research.  For  this  reason,  it  is  crucial  to  investigate  the  factors 
that  have  impact  on  thermal  stresses  in  the  legs,  in  particular  their 
geometry. 

It  is  the  aim  of  this  study  to  investigate  the  possible  effect  of 
thermoelectric  leg  geometries  on  thermo-mechanical  and  power 
generation  performance  of  thermoelectric  devices.  For  this  pur¬ 
pose,  thermoelectric  modules  with  various  leg  geometries 
including  rectangular  prism,  trapezoidal  prism,  cylindrical  and 
octagonal  prism  were  modeled.  Finite-element  analyses  were  car¬ 
ried  out  for  specific  configurations.  Power  generation  through  an 
imposed  steady-state  temperature  gradient  was  considered  as 
operation  mode  of  the  thermoelectric  modules.  Moreover,  modules 
with  segmented  thermoelectric  legs  were  modeled  considering  a 
larger  temperature  gradient.  Temperature  distributions,  magnitude 
of  generated  power  outputs,  conversion  efficiencies,  magnitude 
and  distribution  of  thermal  stresses  were  evaluated. 

2.  Finite  element  modeling 

2.1.  Method 

In  order  to  predict  temperature  distribution,  power  generation 
and  thermal  stresses  with  varying  leg  geometry,  finite-element 
methods  are  used.  Power  generation  and  thermo-mechanical 
behavior  of  a  4-leg  thermoelectric  module  is  considered.  Steady- 
state  temperature  gradients  of  100  °C  and  300  °C  were  applied 
since  most  commonly  used  thermoelectric  devices  are  usually 
designed  to  operate  around  these  temperatures.  Steady-state 
thermal,  thermoelectric,  and  static  structural  modules  of  ANSYS® 
Workbench  (Academic  Research,  Release  14.0)  were  utilized  to 
carry  out  temperature,  power  generation,  and  thermo-mechanical 
analyses,  respectively.  ANSYS  finite  element  software  has  multi¬ 
physics  coupled  simulation  capabilities  to  model  and  analyze 
thermoelectric  effects  with  accurate  solutions  and  is  recognized  by 
the  thermoelectric  device  manufacturers  [30]. 

2.2.  Module  geometries 

Various  thermoelectric  leg  geometries  were  considered.  First,  a 
4-leg  module  with  p-  and  n-type  rectangular  prism  BiTe-based  legs 
with  dimensions  of  1.4  mm  x  1.4  mm  x  3.0  mm  was  modeled.  The 
legs  were  spaced  with  1.0  mm  distance  between  each  other.  This 
module  was  defined  as  the  original  model  since  thermoelectric 
modules  are  typically  made  of  rectangular  prism  legs.  Other  mod¬ 
ules  were  modeled  with  cylindrical  and  trapezoidal-  and 
octagonal-prism  leg  geometries.  Ease  of  fabrication  was  considered 
when  choosing  the  leg  geometries.  Leg  geometries  were  modeled 
with  equal  cross-sectional  areas  in  x-y  plane  and  heights.  Cross- 
sectional  area  at  mid-height  was  considered  for  the  trapezoidal 
prism  model.  Modeled  thermoelectric  modules  with  their  respec¬ 
tive  dimensions  are  shown  in  Fig.  1.  Thicknesses  of  ceramic  sub¬ 
strates,  conductor  pads,  and  solder  strips  are  0.8  mm,  0.15  mm,  and 
0.1  mm  respectively.  Geometries  of  the  copper  pads  and  solder 
strips  were  modified  to  match  leg  geometries.  However,  their 
thicknesses  remained  the  same  for  all  models.  Segmented  leg 
structure  was  considered  for  the  temperature  gradient  of  300  °C. 
Modules  with  rectangular  prism-  and  cylindrical-segmented  legs 
were  modeled.  Segmented  models  share  the  same  dimensions  with 
the  rectangular  prism  and  the  cylindrical  non-segmented  models. 

2.3.  Material  properties 

Thermoelectric  modules  are  usually  composed  of  p-  and  n-type 
thermoelectric  legs  (or  so-called  pins),  ceramic  substrates, 
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octagonal  prism 


cylindrical 


rectangular  prism 


trapezoidal  prism 


Fig.  1.  Models  with  various  leg  geometries  and  dimensions  (in  mm). 


conductor  pads,  Ni-rich  diffusion  barriers,  and  soldering  strips.  BiTe- 
based  alloys  are  the  most  commonly  preferred  leg  materials  for 
relatively  small  temperature  gradient  applications.  In  our  models, 
extruded  p-type  (Bi0.5Sb0.5)2Te3  and  n-type  (Bio.gsSbo.osM- 
Teo.gsSeo.osja  alloys  were  preferred  as  the  thermoelectric  leg  mate¬ 
rials  due  to  their  optimal  performance  for  our  defined  temperature 
range  of  20-120  °C.  In  case  of  large  temperature  gradient  models, 
segmented  legs  were  used  for  optimized  Seebeck  voltage  and 
maximized  power  output.  Initial  segments  of  the  legs  were  modeled 
with  the  extruded  BiTe-based  alloys  and  second  segments  were 
modeled  with  CoSb-based  (Skutterudite)  p-type  Ceo.gFesCoSb^  and 
n-type  Co4SbnGeo.4Teo.6  alloys.  CoSb-based  alloys  were  preferred 
due  to  their  optimal  Seebeck  coefficients  for  high  temperature  ap¬ 
plications.  Temperature-dependent  Seebeck  coefficients  of  both 
extruded  BiTe-  and  CoSb-based  alloys  are  shown  in  Fig.  2.  Other 
material  properties  of  them  are  listed  in  Table  1. 

Since  temperatures  vary  across  the  thermoelectric  device, 
temperature-dependent  properties  of  the  materials-including  the 
Seebeck  coefficient,  thermal  conductivity,  coefficient  of  thermal 
expansion,  and  electrical  resistivity-were  taken  into  account  in  the 
finite-element  models.  The  ANSYS  Workbench  analysis  system 
supports  temperature-dependent  problems  [31].  When  a 
temperature-dependent  property  is  defined,  the  program  includes 
it  in  the  solution  of  element  matrices.  The  corresponding  material 
properties  in  the  element  matrices  are  determined  with  calculation 
of  temperatures  at  the  center  points  of  the  elements  and  linear 


interpolation  of  the  material  property  values  for  these  tempera¬ 
tures.  If  the  calculated  temperatures  are  beyond  the  extreme 
minimum  and  maximum  temperatures,  the  program  determines 
corresponding  material  properties  using  constant-valued 
extrapolations. 


Fig.  2.  Seebeck  coefficients  of  p-type;  (Bi0.2Sb0.8)2Te3  and  Ce0.9Fe3CoSbi2  alloys  and  n- 
type;  (Bio.gsSbo.osMTeo.gsSeo.osh  and  Co4SbiiGeo.4Teo.6  alloys  [18,32-35]. 


U.  Erturun  et  al.  /  Applied  Thermal  Engineering  73  (2014)  126-139 


129 


Table  1 

Properties  of  p-type  (Bio.2Sbo.8)2Te3  and  n-type  (Bio.gsSbo.osMTeo.gsSeo.osh  alloys 
[18,20,32,38,39],  along  with  p-type  Ce0.gFe3CoSbi2  and  n-type  Co4SbnGeo.4Te0.6 
alloys  [27,33,34,40]. 


Property 

p-type 

BiTe-based 

alloy 

n-type 

BiTe-based 

alloy 

p-type 

CoSb-based 

alloy 

n-type 

CoSb-based 

alloy 

d  (kg/m3) 

6858 

7858 

7656 

7539 

a  (1(T6/K) 

16.8 

16.8 

11.5 

11.5 

E  (GPa) 

130 

130 

Radial 

40.3 

42.7 

Axial 

49.7 

51.0 

V 

0.28 

0.28 

0.23 

0.23 

X  (W/m  K) 

1.22-5.24 

1.78-3.50 

1.87-1.77 

4.49-3.20 

(300-500  I<) 

(300-500  I<) 

(300-600  K) 

(300-600  K) 

p  (|i(2  m) 

13.8-25.0 

9.8-14.9 

10.43-12.65 

11.07-13.87 

(300-500  I<) 

(300-500  I<) 

(300-600  K) 

(300-600  K) 

Copper  alloy  was  preferred  as  the  conductor  material  due  to  its 
high  electrical  and  thermal  conductivity.  In  order  to  join  the  legs 
and  conductor  pads  two  different  soldering  alloy  were  used. 
63Sn-37Pb  alloy  was  preferred  for  the  small  temperature  gradient 
models.  For  the  large  temperature  gradient  models, 
92.5Pb-5Sn-2.5Ag  alloy  was  preferred  due  to  its  superior  elec¬ 
trical,  thermal,  and  mechanical  properties.  Moreover  durability  of 
92.5Pb-5Sn-2.5Ag  alloy  is  higher  than  63Sn-37Pb  alloy  at  high 
temperatures.  In  order  to  provide  mechanical  integrity  of  the 
module,  AI2O3  ceramic  substrates  were  preferred.  Material  prop¬ 
erties  of  the  conductors,  soldering  alloy,  and  ceramic  substrates  are 
listed  in  Table  2.  Temperature  dependencies  of  some  material 
properties  including  thermal  and  electrical  conductivities  were 
considered  in  the  simulations.  Thin  Ni-rich  diffusion  barriers  were 
not  modeled  since  their  effect  on  thermo-mechanical  behavior  of 
the  thermoelectric  module  is  negligible  due  to  their  very  small 
thicknesses,  usually  4  pm,  compared  to  thicknesses  of  100  pm  for 
the  solder  and  150  pm  for  the  conductor  parts. 

Research  in  the  literature  revealed  that  thermal  stresses  in  the 
thermoelectric  legs  can  be  significantly  lowered  due  to  the  elastic 
and  plastic  deformation  of  the  conductor  and  soldering  materials 
[18,20,27].  Particularly  soldering  alloys  have  an  important  role  to 
reduce  the  stress  levels  in  the  legs  since  their  yield  strength  is 
lower  than  the  strength  of  leg  material.  Clin  et  al.  [18]  concluded 
that  stresses  in  the  legs  could  be  two  or  three  times  higher  if  the 
plastic  deformation  of  the  soldering  alloy  is  not  considered.  Thus, 
elasto-plastic  behavior  of  the  conductor  and  the  soldering  alloys  is 
considered  in  our  finite-element  models.  Bilinear  kinematic 
hardening  is  preferred  between  the  yield  and  the  ultimate  tensile 
stresses  of  the  Cu  conductor  and  the  63Sn-37Pb  solder.  Elonga¬ 
tion,  yield  stresses,  and  ultimate  tensile  stresses  of  the  Cu 


Table  2 

Properties  of  A1203  (A4,  99.7%),  Cu-based,  63Sn-37Pb,  and  92.5Pb-5Sn-2.5Ag 
alloys  [20,41-48]. 


Property 

ai2o3 

substrate 

Cu-based 

conductor 

63Sn— 37Pb 
solder 

92.5Pb— 5Sn— 2.5Ag 
solder 

d  (kg/m3) 

3970 

8940 

8420 

11300 

a  (10-6/I<) 

4.89-6.68 

16.7-17.6 

24 

29 

(223-600  K) 

(223-583  I<) 

E  (GPa) 

380-366 

119-100 

36-12 

23-15 

(293-600  I<) 

(273-600  I<) 

(273-423  K) 

(300-600  I<) 

V 

0.26 

0.31 

0.40 

0.44 

ay  (MPa) 

70 

20 

(Tuts  (MPa) 

250 

40 

%EL 

69 

40 

X  (W/m  K) 

37.2-17.2 

398-382 

53-46 

35 

(293-600  K) 

(293-600  I<) 

(273-423  K) 

p  (|iQ-m) 

1  x  1018 

0.18-0.38 

0.12-0.23 

0.2 

(300-600  I<) 

(293-473  I<) 

Table  3 

Anand  viscoplasticity  model  constants  of  92.5Pb— 5Sn-2.5Ag  alloy  [37]. 


So  (MPa) 

d/R  (I<) 

Do  (sec  1) 

f  m 

h0  (MPa) 

S  (MPa)  n  a 

33.07 

11024 

1.052(105) 

7  0.241 

1432 

41.63  0.002  1.3 

conductor  and  the  63Sn-37Pb  solder  are  shown  in  Table  2.  Anand 
viscoplasticity  is  a  commonly  used  model  in  finite-element  sim¬ 
ulations  to  simulate  elasto-plastic  behavior  of  92.5Pb-5Sn-2.5Ag 
solders  [36,37].  In  our  large  temperature  gradient  (300  °C)  models, 
Anand  viscoplasticity  was  applied  since  it  is  a  verified  model  for 
deformation  behavior  of  92.5Pb-5Sn-2.5Ag  solders  at  high  tem¬ 
peratures.  The  model  constants  for  92.5Pb-5Sn-2.5Ag  alloy  are 
listed  in  Table  3. 

2.4.  Simulation  parameters  and  meshing 

Thermoelectric  power  generation  was  considered  while 
defining  simulation  parameters  for  the  models.  In  case  of  small 
temperature  gradient  models,  steady  state  temperature  values  of 
120  °C  and  20  °C  were  imposed  on  the  hot  and  cold  sides  of  the 
module,  respectively.  In  case  of  large  temperature  gradient  models, 
imposed  temperature  values  were  320  °C  and  20  °C.  Environmental 
temperature  was  assumed  to  be  20  °C  for  all  models.  Seebeck 
electric  potentials  were  calculated  analytically  using  the  formula 
below 

V  =  S(Th  —  Tc)  (1) 

where  S  is  Seebeck  coefficient,  TH  is  hot  side  and  Tc  is  cold  side 
temperature.  Calculated  Seebeck  electric  potentials  for  the  small 
and  large  temperature  gradient  models  were  45  mV  and  108  mV, 
respectively.  These  electric  potential  values  are  considered  as  high 
electric  potentials  and  applied  to  the  p-type  conductor  pads.  Low 
electric  potentials  are  defined  as  0  mV  for  the  n-type  conductor 
pads.  Air  convection  with  the  coefficient  of  10  W/m2  °C  is  assumed 
and  applied  to  the  outer  surfaces  of  the  modules.  The  modules  are 
mechanically  unconstrained  and  remain  free-standing  during  the 
simulations.  Medium  size  meshing  is  preferred  using  the  automatic 
mesh  feature  of  ANSYS.  Defined  parameters  and  meshed  models 
are  shown  in  Fig.  3(a)  and  (b)  for  the  temperature  gradient  of  100  °C 
and  300  °C  cases,  respectively. 

In  order  to  establish  mesh  convergence,  the  results  of  a  critical 
parameter  (stress,  strain,  or  deflection)  in  a  specific  location  is  be¬ 
ing  plotted  against  a  measure  of  mesh  density.  Simulation  param¬ 
eters  illustrated  in  Fig.  3(a)  are  used.  As  shown  in  Fig.  4,  the 
equivalent  stress  results  of  a  very  critical  point— a  hot  side  corner  of 
a  rectangular  prism  leg— was  plotted  against  the  number  of  nodes 
after  eight  convergence  runs.  This  plot  can  be  used  to  indicate  when 
the  convergence  is  achieved.  The  model  converges  above  46  x  103 
nodes  when  the  stress  value  difference  between  the  runs  remains 
less  than  0.7%.  According  to  these  results,  the  finite-element 
models  used  in  this  study  converge  since  their  mesh  patterns 
contain  at  least  75  x  103  nodes. 

2.5.  Model  verification  method 

In  order  to  validate  our  finite-element  model,  numerical  results 
were  compared  with  experimental  measurements.  Height  differ¬ 
ences  AZ  between  a  corner  and  the  center  point  on  the  hot  side 
surface  of  a  commercially  available  40  x  40  mm2  thermoelectric 
module,  Kryotherm  TB-127-1.4-1.5,  were  experimentally  measured 
for  different  temperature  gradients.  Total  thickness  of  the  module  is 
3.9  mm  and  consists  of  p-type  (Bi0.5Sb0.5)2Te3  and  n-type 
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Fig.  3.  Finite-element  models  and  simulation  parameters  for  the  temperature  gradient  of  (a)  100  °C  and  (b)  300  °C. 


Fig.  4.  Convergence  of  the  predicted  equivalent  stresses  in  a  rectangular  prism  leg  as  a 
function  of  the  total  number  of  nodes  in  its  mesh  pattern. 


(Bi0.95Sb0.05)2(Te0.95Se0.05)3  legs  (total  number  of  256)  with  a 
height  of  1.5  mm  and  a  width  of  1.4  mm.  An  experimental  setup 
consists  of  an  Ambios  XP-1  stylus  profiler,  a  Fluke  52  K/J  digital 
thermometer,  type  I<  thermocouples,  a  laser  thermometer  and  a 
GW  Instek  PSP405  programmable  DC  power  supply  was  used.  The 
profiler  has  maximum  vertical  range  of  400  pm  and  vertical  reso¬ 
lution  of  15  A  at  100  pm.  In  order  to  measure  profiles:  (1)  the 
module  was  heated  by  applying  certain  DC  voltages  and  currents 
(and  thereby  obtaining  the  Peltier  effect),  (2)  temperatures  and 
corresponding  power  inputs  were  recorded,  (3)  thermocouples 
were  removed  and  the  profiler  probe  was  positioned  on  the  cen¬ 
ter-starting  point— of  the  hot  side  surface,  (4)  profiles  were 
measured  by  running  the  stylus  profiler  for  each  temperature 
gradient  while  temperatures  were  monitored  with  the  non-contact 
infrared  thermometer  and  (5)  A Z  between  the  center  and  corner 
points  were  calculated  using  the  profile  data.  Profiler  setup  and 
measured  AZ  are  illustrated  through  Fig.  5(a)— (c). 

To  predict  the  height  differences,  thermo-mechanical  finite 
element  analyses  of  a  quarter  model— since  the  module  is  quasi- 
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Fig.  5.  Profile  measurements;  (a)  inside  of  the  stylus  profiler,  (b)  direction  of  the  probe  and  (c)  measured  height  difference  (AZ). 


symmetric— were  carried  out.  Our  4-leg  finite-element  model 
showed  in  Fig.  3(a)  was  modified  to  a  quasi-symmetric  64-leg 
model  has  matching  materials  and  dimensions  of  the  quarter  of 
the  Kryotherm  device.  Thermal  and  mechanical  material  properties 
of  p-  and  n-type  BiTe-based  legs,  92.5Pb-5Sn-2.5Ag  solder, 
Cu-based  conductor  and  AI2O3  substrate,  listed  in  Tables  1  and  2, 
were  used. 

3.  Results  and  discussions 

3.1.  Temperature  analysis 

Temperature  distributions  of  the  thermoelectric  modules  were 
calculated  through  steady-state  thermal  analysis.  Imposed  tem¬ 
perature  gradient  causes  heat  flow  through  entire  module.  For 
temperature  gradient  of  100  °C,  temperature  distributions  of  a 
thermoelectric  element  (coupled  p-  and  n-type  legs)  with 
conductor  pads  and  soldering  strips  of  each  model  are  presented 
in  Fig.  6(a).  Other  parts  of  the  thermoelectric  modules  are  hidden 
in  the  figure  for  only  visualization  purposes.  In  addition,  temper¬ 
ature  distributions  along  the  defined  paths  in  the  z-axis  shown  in 
Fig.  6(a)  are  plotted  and  presented  in  Fig.  6(b).  Temperature  decays 
significantly  along  the  z-axis  between  the  hot  and  cold  side  while 
it  remains  almost  equal  along  the  x-  and  y-axis  for  the  all  models 
as  shown  in  Fig.  6(b).  All  models  except  the  trapezoidal  prism 
presented  fairly  similar  and  linear  temperature  profile  changes 
along  the  z-axis  path;  trapezoidal  prism  legs  presented  slightly 
different  temperature  decays  than  the  other  legs.  This  result  is  a 


consequence  of  various  cross-sectional  areas  of  the  trapezoidal 
prism  legs. 

Temperature  distributions  of  the  segmented  models  for  the 
temperature  gradient  of  300  °C  are  shown  in  Fig.  7(a)  and  (b). 
Rectangular  prism-  and  cylindrical-segmented  legs  presented  a 
fairly  similar  temperature  distribution  trend.  Temperature  decays 
along  the  z-axis  were  not  linear  in  this  case  as  shown  in  Fig.  7(b). 
Noticeable  changes  in  the  slope  of  the  profiles  can  be  clearly  seen 
above  the  mid-height  of  the  legs.  This  non-linearity  occurs  due  to 
the  different  thermal  conductivities  of  BiTe-  and  CoSb-based  alloys. 
Overall,  thermal  analysis  results  show  that  temperature  distribu¬ 
tions  were  slightly  affected  by  changing  thermoelectric  leg 
geometries. 

3.2.  Power  generation  analysis 

In  order  to  evaluate  the  effect  of  thermoelectric  leg  geometries 
on  the  power  generation  performance  of  thermoelectric  modules, 
finite-element  thermoelectric  analyses  were  performed  on  rect¬ 
angular-,  trapezoidal-,  octagonal-prism,  and  cylindrical  models  for 
the  temperature  gradient  of  100  °C.  Applied  temperature  gradients 
generate  electric  current  flow  through  p-  and  n-type  legs  in  series. 
Total  current  densities  shown  as  contours  in  Fig.  8  present  the 
current  flows  in  various  leg  geometries.  Electric  current  and 
absorbed  heat  magnitudes  were  obtained  as  finite-element  anal¬ 
ysis  results.  Generated  power  outputs  and  conversion  efficiencies 
were  analytically  calculated  applying  the  numerical  results,  the 
electric  current  and  absorbed  heat,  in  the  equations  below 
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(a) 


Fig.  6.  Temperature  distributions  in  the  legs  (a)  as  nephogram  and  (b)  along  the  paths  for  the  temperature  gradient  of  100  °C. 
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Fig.  7.  Temperature  distribution  in  the  segmented  legs  (a)  as  nephogram  and  (b)  along  the  paths  for  the  temperature  gradient  of  300  °C. 


P=,n[s(th-Tc)i  -I2r]  (2) 

R  =  A(Pn  +  pp)  (3) 


where  P  is  output  power,  N  is  number  of  thermoelectric  elements  in 
the  module  (in  our  models;  N  =  2),  S  is  Seebeck  coefficient,  Th  is  hot 
side  temperature,  Tc  is  cold  side  temperature,  I  is  electric  current,  R 
is  combined  internal  electrical  resistance  for  the  p-  and  n-type 
thermoelectric  legs,  A  and  L  are  cross-sectional  area  and  length  of 
the  legs,  pp  and  pn  are  electrical  volume  resistivity  for  the  p-  and  n- 
type  legs  respectively,  7]  is  conversion  efficiency  and  Qh  is  absorbed 
heat.  Joule  heating  I2R  due  to  internal  electrical  resistance  of 


thermoelectric  materials  [49]  was  considered  in  the  power  calcu¬ 
lation.  Generated  electric  current  I  and  absorbed  heat  Qh  as  finite- 
element  analysis  results  and  calculated  conversion  efficiencies  7]  are 
shown  in  Fig.  9.  Calculated  power  outputs  for  rectangular  prism, 
trapezoidal  prism,  cylindrical,  and  octagonal  prism  models  were 
18.31  mW,  17.64  mW,  18.85  mW,  and  18.54  mW,  respectively. 
Although  the  electric  current  flows  were  slightly  affected  by 
geometrical  differences,  the  changes  in  the  power  outputs 
remained  very  small.  This  result  was  expected  since  the  parameters 
including  Seebeck  coefficients,  temperature  gradients,  and  heights 
of  the  thermoelectric  legs  were  the  same  for  all  models.  Disparities 
in  power  output  and  efficiency  values  of  the  trapezoidal  prism 
model  were  consequence  of  its  various  cross-sectional  leg  areas. 

In  order  to  calculate  electrical  current  and  absorbed  heat  for  the 
temperature  gradient  of  300  °C,  numerical  thermoelectric  analyses 
were  performed  on  rectangular  prism-  and  cylindrical-segmented 
models.  Power  outputs  and  conversion  efficiencies  were  calculated 
using  finite-element  analysis  results  (/  and  Qh)  and  equations 
(2)— (4).  Generated  electric  current  /,  absorbed  heat  Qh  and 
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Fig.  8.  Total  current  densities  for  the  temperature  gradients  of  (a)  100  °C  and  (b)  300  °C. 


calculated  conversion  efficiencies  rj  are  shown  in  Fig.  9  for  the 
temperature  gradient  of  300  °C.  Calculated  power  outputs  for  rect¬ 
angular  prism-segmented  and  cylindrical-segmented  models  were 
143.6  mW  and  143.7  mW  respectively.  It  can  be  seen  that,  power 
outputs  were  almost  8  times  greater  and  conversion  efficiencies 
were  almost  2  times  greater  than  those  of  the  temperature  gradient 


of  100  °C  case.  This  is  as  a  result  of  higher  Seebeck  coefficients  in 
segmented  models  as  well  as  larger  temperature  gradients.  Perfor¬ 
mance  difference  between  rectangular  prism-segmented  and 
cylindrical-segmented  models  remained  negligible.  It  is  noted  that, 
optimization  of  the  power  output  performance  or  conversion  effi¬ 
ciency  is  not  the  aim  of  this  study. 
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Fig.  9.  Predicted  electric  current,  absorbed  heat,  and  calculated  conversion  efficiencies. 


3.3.  Thermo-mechanical  analysis 

In  order  to  evaluate  the  impact  of  various  thermoelectric  leg 
geometries  on  thermo-mechanical  performance,  thermal  stresses 
developed  in  the  legs  due  to  temperature  gradients  and  thermal 
expansion  mismatches  were  analyzed  for  steady-state  condi¬ 
tions.  Coupling  temperature  and  static-structural  analyses 
compose  thermo-mechanical  analysis.  Temperature  distributions 
obtained  from  the  thermal  analysis  were  applied  as  load  for 
static-structural  analyses.  Thermal  stresses  were  developed  in 
the  legs  due  to  expansion  of  the  hot  side  ceramic  substrate  while 
the  cold  side  obstructs  and  thermal  expansion  mismatches  be¬ 
tween  the  materials  constituting  the  thermoelectric  module. 
Thermal  stress  levels  were  calculated  in  terms  of  equivalent 
stresses  ov  (so-called  von  Mises  stresses).  Equivalent  stress  is 
preferred  since  it  is  an  average  scalar  stress  value  computed  from 
the  second  deviatory  stress  components  o-y,  cauchy  or  so-called 
true  stress  tensors  as  formulated  below: 


(°Tl  -0-22)2  +  (cT22  — (T33)2  +  (<733  _(Tll)2  +  6((7|2  +  (T^3+(7^1) 

2 

(5) 

Distribution  and  maximum  values  of  calculated  equivalent 
stresses  in  the  legs  for  rectangular  prism,  trapezoidal  prism,  cy¬ 
lindrical,  and  octagonal  prism  models  are  shown  in  Fig.  10(a)  for 
temperature  gradient  of  100  °C.  As  expected,  the  maximum  stress 
zones  were  developed  in  the  legs  on  the  hot  side  of  the  module. 
This  is  due  to  higher  thermal  expansions  and  larger  CTE  differences 
between  the  materials  on  the  hot  side.  Maximum  stresses  occur  on 
the  hot  side  corner  of  the  legs  (for  cylindrical  models  on  the  edges) 
and  their  values  for  rectangular  prism,  trapezoidal  prism,  cylin¬ 
drical,  and  octagonal  prism  models  were  49.9  MPa,  41.5  MPa, 
43.3  MPa,  and  47.8  MPa,  respectively.  In  order  to  evaluate  variation 
of  the  stress  values  from  one  leg  to  the  other  in  different  directions, 
stress  paths  taking  the  maximum  stress  node  as  starting  point  were 
constructed.  Equivalent  stress  distributions  along  the  paths  in  the 
x-,  y-,  and  z-axis  are  illustrated  in  Fig.  10(b),  (c)  and  (d).  We  observe 
that  the  stress  values  were  smaller  in  the  middle  section  and  larger 
on  the  edge  of  the  legs.  Along  the  path  of  x-axis,  stresses  distributed 
more  uniformly  in  the  trapezoidal  prism  legs  compared  to  the 
others.  This  can  be  attributed  to  the  longer  hot  side  edges  (1.8  mm 
compared  to  1.4  mm  of  the  rectangular  prism  leg)  of  the  trapezoidal 
prism  legs.  Along  the  path  of  z-axis,  thermal  stress  levels  increase 
closer  to  the  mid-height  sector  along  with  the  hot  and  cold  side 
surfaces.  Stress  values  along  the  paths  ofy-  and  z-axis  are  smaller  in 


the  cylindrical  and  the  octagonal  prism  legs  compared  to  the 
others.  Maximum  principal  stresses  in  the  legs  were  also  calcu¬ 
lated:  40.5  MPa,  42.5  MPa,  33.3  MPa,  and  35.7  MPa  for  rectangular 
prism,  trapezoidal  prism,  cylindrical,  and  octagonal  prism  models, 
respectively.  These  results  show  that  changing  leg  geometry  can 
significantly  reduce  and  redistribute  thermal  stresses  in  the  ther¬ 
moelectric  legs.  Cylindrical  or  trapezoidal  prism  legs  can  be 
preferred  over  the  other  leg  geometries  due  to  their  small  stress 
levels  and  ease  of  fabrication  (particularly  with  respect  to  the  dicing 
process)  compared  to  the  octagonal  prism  models.  To  fabricate  a 
typical  thermoelectric  device,  thermoelectric  legs  are  diced  in  the 
desired  dimensions  from  the  plated  p-  and  n-type  semiconductor 
wafers  using  regular  cutting  saws  or  Electrostatic  Discharge 
Machining  (EDM).  Copper  conductor  pads  and  AI2O3  ceramic  sub¬ 
strates  are  also  diced  to  the  proposed  dimensions. 

Predicted  maximum  plastic  strains  on  the  conductor  pads  were 
0.33%,  0.32%,  0.24%,  and  0.25%  for  rectangular  prism,  trapezoidal 
prism,  cylindrical,  and  octagonal  prism  models,  respectively. 
Maximum  plastic  strains  on  the  soldering  strips  were  0.57%,  0.67%, 
0.65%,  and  0.74%  for  rectangular  prism,  trapezoidal  prism,  cylin¬ 
drical,  and  octagonal  prism  models,  respectively.  These  values  are 
much  smaller  than  the  plastic  extensibilities  of  69%  for  the  copper 
conductor  and  40%  for  63Sn— 37Pb  soldering  alloys  shown  in 
Table  2.  It  is  clear  that  both  conductor  pads  and  soldering  strips 
were  plastically  deformed,  but  they  were  still  in  the  safe  zone  and 
will  not  fail  for  the  defined  conditions.  In  our  previous  thermo¬ 
mechanical  analysis  of  thermoelectric  modules  [50],  only  elastic 
behavior  of  the  soldering  and  conductor  alloys  (zero  plastic  defor¬ 
mation)  was  taken  into  account.  The  maximum  equivalent  stress  of 
129  MPa  was  calculated  in  the  rectangular  prism  thermoelectric 
legs  for  the  temperature  gradient  of  138  °C.  In  the  present  study,  the 
maximum  equivalent  stress  of  49.9  MPa  was  calculated  in  the 
rectangular  prism  legs  for  the  temperature  gradient  of  100  °C.  This 
decrease  in  the  stress  levels  occurred  due  to  the  elastic  and  plastic 
deformations  of  the  conductor  and  soldering  alloys.  Elasto-plastic 
behavior  of  the  conductor  and  soldering  alloys  is  more  appro¬ 
priate  to  represent  the  stress  levels  compare  to  the  linear  elastic 
behavior. 

In  case  of  temperature  gradient  of  300  °C,  developed  thermal 
stresses  in  the  legs  of  the  rectangular  prism-  and  cylindrical- 
segmented  models  were  analyzed.  Levels  and  distribution  of 
calculated  equivalent  stresses  are  shown  in  Fig.  11(a).  Maximum 
thermal  stress  magnitudes  were  smaller  in  the  cylindrical 
segmented  legs  compare  to  the  stresses  in  the  rectangular  prism 
segmented  legs.  In  the  rectangular  prism  segmented  model,  a 
calculated  maximum  equivalent  stress  of  94.2  MPa  was  occured 
between  the  CoSb-  and  BiTe-based  segments.  In  the  cylindrical 
segmented  model,  maximum  equivalent  stress  of  85.7  MPa  was 
occured  on  the  hot  side  edge  of  the  CoSb-based  segment.  Distri¬ 
bution  of  the  stresses  is  illustrated  by  using  the  paths  shown  in 
Fig.  11(a)  and  the  profiles  along  these  paths  in  Fig.  11(b),  (c),  and 
(d).  Stress  evolution  of  the  segmented  legs  displays  distinctive 
trends  compared  to  the  non-segmented  models.  Significant 
changes  in  the  stress  magnitudes  and  distributions  were  observed 
in  the  joints  of  the  segments  as  shown  along  the  path  of  z-axis  in 
Fig.  11(d).  In  the  segmented  models  high  stress  levels  occur  be¬ 
tween  the  segments  due  to  the  thermal  expansion  differences 
between  BiTe-  and  CoSb-based  alloys.  It  is  noted  that  different 
soldering  alloy  and  plasticity  models  were  used  in  the  segmented 
models.  This  made  impact  on  stress  trends  on  the  hot  side  of  the 
legs  as  shown  along  the  path  ytop  in  Fig.  11(c).  Maximum  plastic 
strains  of  1%  and  3.4%  were  developed  in  the  conductor  pads  and 
soldering  strips  of  the  rectangular  prism-  segmented  model.  They 
were  0.93%  and  0.64%,  respectively  for  the  cylindrical-segmented 
model.  Overall,  increases  in  thermal  stresses  magnitude  (due  to 
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Fig.  10.  Magnitude  and  distribution  of  maximum  equivalent  stresses  in  the  legs;  (a)  as  nephogram  and  along  the  (b)  x-,  (c)  y-,  and  (d)  z-axis  paths  for  the  temperature  gradient  of 
100  °C. 
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Fig.  11.  Maximum  equivalent  stresses  in  the  segmented  legs;  (a)  as  nephogram  and  along  the  (b)  ymid-  (c)  ytoP-.  and  (d)  z-axis  paths  for  the  temperature  gradient  of  300  °C. 


larger  temperature  gradients)  and  distinctive  stress  distributions 
were  observed  in  the  segmented  legs. 

3.4.  Verification  results 

Experimental  profile  measurements  and  numerical  simulations 
were  performed  using  the  temperatures,  shown  in  Table  4.  Height 
differences  AZ  between  the  center  point  and  the  corner  of  the  hot 
side  of  the  module  were  experimentally  measured  for  the  defined 
temperatures.  Corresponding  AZ  values  were  numerically 


Table  4 

Measured  hot  and  cold  side  temperatures  of  the  Kryotherm  thermoelectric  device. 


# 

Th  (°C) 

Tc(°C) 

AT(°C) 

1 

49 

32 

17 

2 

61 

40 

21 

3 

70 

44 

26 

4 

80 

48 

32 

5 

90 

63 

27 

6 

100 

73 

27 

7 

110 

80 

30 

8 

122 

91 

31 
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Fig.  12.  Results  show  (a)  numerical  simulation  of  deflection  and  (b)  comparison  of  experimental  and  numerical  AZ  values. 


predicted  using  the  64-leg  quasi-symmetric  model  for  the  same 
temperatures  and  the  boundary  conditions.  Results  of  ana¬ 
lysis-temperature  distributions  and  mechanical  deforma¬ 
tions— are  shown  in  Fig.  12.  Moreover,  our  numerical  results  were 
also  compared  to  results  of  a  study  in  the  literature.  Turenne  et  al. 
[20]  modeled  a  Kryotherm  TB-127-1.4-1.50  thermoelectric  device 
and  predicted  a  AZ  of  73  pm  between  the  center  and  the  external 
corner  of  the  hot  side  for  the  temperature  gradient  of  100  °C.  They 
experimentally  measured  a  corresponding  AZ  of  63  pm  using  a 
laser  profiler.  In  our  finite-element  model,  a  AZ  of  73.74  pm  was 
predicted  for  the  temperature  gradient  of  100  °C  under  the  same 
boundary  conditions  with  their  simulations.  These  results  validate 
our  finite-element  model  for  the  thermo-mechanical  analysis. 


4.  Conclusion 

Effect  of  thermoelectric  leg  geometries  on  power  generation 
and  thermo-mechanical  performance  of  thermoelectric  devices 
was  evaluated  utilizing  numerical  methods.  Rectangular  prism, 
trapezoidal  prism,  cylindrical,  and  octagonal  prism  leg  geome¬ 
tries  and  steady-state  temperature  gradients  of  100  °C  and 
300  °C  were  considered.  Thermoelectric  modules  with  BiTe- 
based  non-segmented  and  BiTe/CoSb-based  segmented  legs 
were  modeled.  Temperature-dependent  material  properties  and 
stress  relaxation  mechanisms  due  to  elasto-plastic  behavior  were 
considered  in  the  finite-element  models.  The  following  conclu¬ 
sions  were  reached:  (1)  Changes  in  temperature  distributions, 
power  outputs,  and  conversion  efficiencies  were  limited  for 
various  leg  geometries.  (2)  Elastic  and  plastic  deformations  of  the 
soldering  and  conductor  alloys  helped  to  decrease  the  stress 
levels  in  the  legs.  (3)  Changes  in  thermal  stress  distributions  and 
magnitude  of  maximum  equivalent  stresses  were  significant  for 
various  leg  geometries.  Maximum  equivalent  stresses  in  the 
rectangular  prism  and  cylindrical  legs  were  49.9  MPa  and 
43.3  MPa,  respectively  for  the  temperature  gradient  of  100  °C.  (4) 
Maximum  equivalent  stresses  were  94.2  MPa  and  85.7  MPa  in  the 
rectangular  prism-segmented  and  cylindrical-segmented  legs  for 
the  temperature  gradient  of  300  °C.  In  the  segmented  models, 
high  stresses  were  occurred  also  between  the  segments  along 
with  the  hot  sides  of  the  legs.  (5)  It  is  clear  that  in  both  non- 
segmented  and  segmented  models,  leg  geometries  have  a 
significant  effect  on  the  thermo-mechanical  performance  of 
thermoelectric  devices.  In  order  to  reduce  thermal  stresses 
developed  in  thermoelectric  devices,  varying  thermoelectric  leg 
geometry  can  be  considered  as  a  factor. 


Nomenclature 

S  Seebeck  coefficient,  V/I< 

V  Seebeck  voltage,  V 

Qh  hot  side  heat  absorbed,  W 

P  useful  output  power,  W 

N  number  of  thermoelectric  elements,  — 

R  electrical  resistance  for  a  thermoelectric  element,  pQ 

L  length  of  a  thermoelectric  leg,  m 

A  cross-sectional  area  of  a  thermoelectric  leg,  m2 

I  electric  current,  A 

Th  hot  side  temperature,  I< 

Tc  cold  side  temperature,  K 

d  density,  kg/m3 

a  coefficient  of  thermal  expansion,  K-1 

X  thermal  conductivity,  W/m  I< 

E  modulus  of  elasticity,  GPa 

v  poison's  ratio,  - 

Gy  yield  strength,  MPa 

(Juts  ultimate  tensile  strength,  MPa 

%EL  percent  elongation  at  yield,  - 

p  electrical  volume  resistivity,  pQ  •  m 

So  initial  value  of  deformation  resistance,  MPa 

Q.  activation  energy,  ev 

R  ideal  gas  constant,  J/kg  K 

Do  pre-exponential  factor,  sec  1 

?  multiplier  of  stress,  - 

m  strain  rate  sensitivity  of  stress,  - 

ho  hardening/softening  constant,  MPa 

S  coefficient  for  deformation  resistance  saturation,  MPa 

n  strain  rate  sensitivity  of  deformation  resistance,  - 

a  strain  rate  sensitivity  of  hardening  or  softening,  - 

r i  efficiency,  - 

Gy  stress  tensor  components,  MPa 

Gy  equivalent  (von-Mises)  stress,  MPa 

AZ  height  difference,  pm 
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